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COLLIDING ELECTROMAGNETIC WAVES: 
EXACT SOLUTIONS 

It is the purpose of this chapter and the following one to review the exact 
solutions which describe colliding plane electromagnetic waves or com- 
binations of gravitational and electromagnetic waves that are presently 
known. Since these space-times have two Killing vectors and an Ernst 
method can be used, there are a whole range of generating techniques 
that can be applied to obtain new classes of solutions. It is therefore 
not possible to review exact solutions without also reviewing the gener- 
ating techniques that have been applied to obtain them. The techniques 
described may clearly be used to obtain further solutions. This chapter 
reviews classes of non-diagonal solutions. The particular case of diagonal 
solutions is considered in Chapter 18. 

17.1 A technique of Chandrasekhar and Xanthopoulos 

Having developed the approach described in the previous chapter, Chan- 
drasekhar and Xanthopoulos (1985 a) then used it to obtain an electro- 
magnetic generalization of the Nutku-Halil solution. 

If the Khan-Penrose and Nutku-Halil solutions are regarded as the 
analogues of the Schwarzschild and Kerr solutions respectively, then it is 
reasonable to suppose that the Kerr-Newman solution may also have a 
colliding wave analogue. Accordingly, it is appropriate 1 to consider the 
special class of solutions for which 

H = Q(Z+1) (17.1) 

where Q is a constant that can always be chosen to be real. In this case, 
(16.20) and (16.21) both reduce to the same equation: 

\ ((1 - 2Q 2 )(Z + Z)- 2Q 2 (ZZ + 1)) (((1 - t 2 )Z t ) t - ((1 - z 2 )Z z ) z ) 
= (1 - 2Q 2 (Z + 1)) ((1 - t 2 )Z 2 - (1 - z 2 )Z 2 ) . 

(17.2) 

1 See Chandrasekhar (1983) section 110. 
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It is also convenient to consider the restriction (17.1) in terms of the asso- 
ciated Ernst potentials E and r] given by (16.23,24). It can immediately 
be seen that, in this case, 

V = 2Q (17.3) 
so that equation (16.25) is automatically satisfied, and (16.25) becomes 



(1 - 4Q 2 - EE) (((1 - t 2 )E t ) t - ((1 - z 2 )E z ) 



= -2E((l-t 2 )E 2 -(l-z 2 )E 2 ) 



(17.4) 



which is virtually identical to (11.19). 

In this case, however, it is important to notice that the functions Z 
and E do not contain the metric functions as in Chapter 11. Rather, 
through (16.19), they contain the potentials described in the previous 
chapter. These are now given in terms of E by 

(1 - AO 2 - EE) (E-E) 2Q 

$ = —i - H = — — (17 5) 

(1-E)(1-E)' \l-E){l-EY 1-E' { } 

Chandrasekhar and Xanthopoulos (1985a) have made use of the im- 
portant result that, if E Q is a solution of the vacuum Ernst equation 
(11.14), then 

E = aE Q (17.6) 
is a solution of the equation (17.4) for an Einstein-Maxwell field, where 



a = - 4Q 2 . (17.7) 

In view of the possible transformation (12.9), there is no loss of generality 
in assuming that a is real and that < Q < |. This technique can 
therefore be characterized by the Ernst potentials 



E = aE , 77 = VI - a 2 . (17.8) 

This result can also be stated in the equivalent form that, if Z Q is 
a solution of the vacuum Ernst equation (11.8), then a new solution of 
the Ernst equations (16.20) and (16.21) for an Einstein-Maxwell field is 
given by 

(l-„) + (l + < ,)3, CT(Z + 1) 

(l + a) + (l-a)Z„ (1 + a) + (1 - a)Z 1 ' 
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It follows immediately from (17.9) that the potentials for the new 
solution are given in terms of the vacuum potentials Z Q = ty Q + i$ G by 

2a 2 . „ - , 4a 2 

U il 17.10) 

2ia , „ - s 4a 



where 



O 2 = |(l + a) + (l-a)Z | 2 



= (1 + a) 2 + 2(1 - a 2 )v|/ + (1 - af^l + $ 2 ). (17 ' n) 

Clearly, it is possible to use this method to obtain new solutions of the 
Einstein-Maxwell equations starting with any vacuum solution describing 
colliding gravitational waves. The new solutions contain an arbitrary 
parameter a. Moreover, the new solutions must reduce to the original 
vacuum solutions in the limit as a — > 1. 

By substituting (17.8) into (16.28), it may be observed that the 
boundary conditions for colliding gravitational and electromagnetic waves 
of this type are automatically satisfied if they were satisfied for the original 
vacuum solution given by E Q or Z Q , and with the same values for k± and 
/c 2 . This implies that the expressions for f(u) and g(v) are unchanged, 
and thus 

U = U Q . (17.12) 

Using the approach described above with (17.1) and (17.6), the met- 
ric function \ can immediately be obtained from (16.14). It is then neces- 
sary to integrate the equations (16.15) to obtain uj. Finally it is necessary 
to integrate (7.9) to determine the remaining metric function M. 

It may be noticed, however, that the transformation (17.9) is con- 
tained in the 'Elders-Harrison' transformation for a particular choice of 
parameters (Harrison, 1968). In this case, it is possible to find a general 
transformation for M, so that an integration of (7.9) is not essential. But 
there does not appear to be any way of avoiding the integration of (16.15) 
to obtain uj. 



17.2 Two particular examples 

In this section two results are presented that have been obtained using the 
above technique. The first of these is a generalization of the Nutku-Halil 
solution and the second a generalization of the vacuum Chandrasekhar- 
Xanthopoulos solution. 
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Chandrasekhar and Xanthopoulos (1985a) have obtained a new col- 
liding wave solution using the above method and taking the initial vacuum 
solution to be the Nutku-Halil solution described in Section 13.1. In this 
case, we have Z Q = \I/ + i$ Q where 



\ 1 — p 2 t 2 — q 2 z< 
2q (l-z 2 )(l-pt) 
p (1 — p 2 t 2 — q 2 z 2 ) 



(17.13) 



These expressions may immediately be deduced from equations (13.11) 
and (13.17) after noting that the vacuum Chandrasekhar-Xanthopoulos 
and the Nutku-Halil solutions are related simply by interchanging the 
expressions for \I/ and $ with those for \ and 10. In this case, it may also 
be noticed that a constant term in $ Q has been altered. 

The metric functions for this solution can now be obtained from 
(16.14) and by integrating equations (16.15) and (7.9). The integration 
process is described in great detail by Chandrasekhar and Xanthopoulos 
(1985 a). The resulting expressions, which are the metric functions of the 
line element (11.4), may be written in the form 

e -u = ^i-p^l-z 2 = l-u 2 -v 2 

_ O 2 (l-p 2 t 2 -q 2 z 2 ) 
X ~ 4a 2 ((1-pt) 2 + q 2 z 2 ) 

__ qz (4ap 2 + (1 - a) 2 [l + q 2 z 2 + (3 - z 2 )(l - pt) 2 ] ) (17-14) 
U ~ 2a 2 p 2 ((1 -pt) 2 + q 2 z 2 ) 

- M _ & (l- P 2 t 2 -q 2 z 2 ) 

4a 2 (l_t2)l/4(!_ z 2)l/4 

where Q is given in terms of the expressions (17.13) by (17.11). 

The structure of this solution has been described in detail by Chan- 
drasekhar and Xanthopoulos (1985a), with particular emphasis on the 
colinear case when q = and p = 1. As usual, there is a strong curvature 
singularity in region IV on the hypersurface u 2 +v 2 = 1, and there are fold 
singularities in regions II and III when u = 1 and v = 1. The approaching 
waves are a combination of step gravitational and electromagnetic com- 
ponents, and the wave fronts also include impulsive gravitational waves. 

A second example of the application of the technique described in 
Section 17.1 is contained in part I of another substantial paper of Chan- 
drasekhar and Xanthopoulos (19876). This is obtained by taking the 
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initial solution to be the vacuum Chandrasekhar-Xanthopoulos (19866) 
solution described in Section 13.3. 

In this second example, the expression for Z Q is taken to be (13.10) 
so that 

* = 1 ~ P 2 t 2 ~ 2qz 

° (l-pt) 2 + q 2 z 2 ' ° (1-pt) 2 + q 2 z 2 ' K ' 

Equations (16.15) may now be solved as in the previous example, and the 
resulting solution may be expressed in terms of the metric functions of 
the line element (11.4) as 

e ~ u = y/i-t 2 Vl-z 2 = l-u 2 -v 2 

VT^VT^ ( (1 - apt) 2 + a 2 q 2 z 2 



a 2 \ 1 — p 2 t 2 — q 2 z 2 

q(l-z 2 ) ( l + a 2 -2apt \ ( 17 - 16 ) 

U = 2 1 272 2T2 + COnSt 

a z p \1 — P t — q z £ J 
e _ M = (1 - apt) 2 + a 2 q 2 z 2 
a 2 \Jl — u 2 \J\ — v 2 

Chandrasekhar and Xanthopoulos (19876) have analysed this solu- 
tion and have shown that it is of algebraic type D. They have further 
shown that it is in fact the Kerr-Newman space-time interior to its ergo- 
sphere. In this sense, it can be considered to be the natural generalization 
of the vacuum Chandrasekhar-Xanthopoulos (19866) solution, which is a 
Kerr space-time interior to the ergo-sphere. 

The singularity structure of this solution is the same as that of 
the vacuum Chandrasekhar-Xanthopoulos solution as described in Sec- 
tion 13.3. The space-like surface on which u 2 + v 2 = 1 is an unstable 
quasiregular singularity rather than a curvature singularity. 

In both of the solutions presented in this section there are impulsive 
gravitational wave components along the wavefronts u = and v = 0. 
Hence these solutions do not satisfy the conditions for Tipler's theo- 
rem 8.1. 

The global structure of this solution has been further considered by 
Hay ward (1989 a) as an extension of his analysis of the degenerate Ferrari- 
Ibanez solution. 



17.3 Another type D solution 



It has already been pointed out in Section 16.3 that the Bell-Szekeres 
solution can be obtained from the vacuum Chandrasekhar-Xanthopoulos 
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solution using the simple generation technique described in Section 16.2. 
Chandrasekhar and Xanthopoulos (1987a) have subsequently obtained a 
generalization of the Bell-Szekeres solution by applying an Ehlers trans- 
formation to the Ernst potential that is used to generate the Bell-Szekeres 
solution. The resulting space-time describing the interaction region of the 
new solution is of algebraic type D. 

It may first be recalled that, if Z OQ is a solution of the vacuum Ernst 
equation (11.15), then another solution is given by 



= f^oo (VI. VI) 

where the parameter of this Ehlers transformation has been altered from 
that of (12.13) by putting c = —2/3. In terms of the associated Ernst 
function, the transformation (12.14) may be restated in the form that, if 
E OQ is a solution of equation (11.14), then a new vacuum solution can be 
obtained using the new potential 

E QO - 0(1 + E QO ) 

E ° = TT#W ( 7 8) 

Now, using the simple generation technique described in Section 16.2, 
an electromagnetic solution of the equations (16.20) and (16.21) is given 
by 

Z = 1 and H = E Q (17.19) 

where E Q is a solution of the vacuum Ernst equation (11.14). The gen- 
eralized non-diagonal Bell-Szekeres solution is expressed in this way in 
Section 16.3 with E Q = pt + iqz. A further generalization of the Bell- 
Szekeres solution may now be obtained by applying the Ehlers transfor- 
mation (17.18) to this expression for E Q . Accordingly, we now consider 
an electromagnetic solution given by the potentials 

Z = 1, H = E oo~^ + E ) 

where 

E 00 =pt + iqz, and p 2 + q 2 = 1. (17.21) 



From these expressions it follows that 

\-p 2 t 2 -q 2 z 2 
(l-[3qz) 2 + [3 2 (l+pt)^ 



1 _ B 2 -2 _ 2 2 

r , V J, ^ $ = - ( 17 - 22 ) 
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The expression for x is now determined by (16.14), and the remaining 
metric functions for this solution can be obtained by integrating equations 
(16.15) and (7.9). The details of the integration process are described by 
Chandrasekhar and Xanthopoulos (1987a), and the resulting expressions 
may be written in the form 

e~ u = Vl-tVl-^ 2 = 1 - u 2 - v 2 

VT^VT^ 2 ((1 - Pqz) 2 +/3 2 (1+ pt) 2 ) 
* (l-p 2 t 2 -q 2 z 2 ) (17 23) 

_ p(l - t 2 ) (1 - 2(3qz + 2(3 2 (1 +p)) + 2/3 V(l -t)(l- z 2 ) 

g(l — p 2 t 2 — q 2 z 2 ) 

where the arbitrary constant in the integral for uj has been chosen such 
that u = when t = 1. 

It may be noticed that this solution reduces to the conformally flat 
Bell-Szekeres solution in the limit as (3 — > 0. It may also be observed that 
it does not contain the special case of a diagonal solution in which u = 0, 
except in this Bell-Szekeres limit. 

The properties of this solution have been described in detail by Chan- 
drasekhar and Xanthopoulos (1987a). They have evaluated the compo- 
nents of the curvature tensor and shown that the space-time in the in- 
teraction region is of algebraic type D, with the repeated principal null 
directions of the Weyl tensor being aligned with the principal null direc- 
tions of the electromagnetic field. The approaching waves in regions II 
and III are a combination of gravitational and electromagnetic waves, and 
the usual fold singularities in these regions still occur. 

One particularly significant feature of this solution is that the surface 
in region IV on which t = 1 does not correspond to a curvature singularity. 
This surface, on which the two waves mutually focus each other, is a 
quasiregular singularity on which the components of the curvature tensor 
remain bounded. In this case, there are impulsive gravitational wave 
components along the wavefronts u = and v = 0, and hence the solution 
does not satisfy the conditions for Tipler's theorem 8.1. In this feature 
it is similar to the Bell-Szekeres solution, the degenerate Ferrari-Ibahez 
solution of Section 10.5 and the vacuum Chandrasekhar-Xanthopoulos 
solution described in Section 13.3. 

The solution given in this section has also been independently ob- 
tained and described by Halilsoy (1988a). Halilsoy has also used the 
simplification (17.19) as described in Section 16.2, but with the initial 
solution E being that obtained from a colinear solution using the trans- 
formation (12.44) and (12.45) with a = 1. However, as is clear from 
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(12.46), this can also be considered as an Ehlers transformation (12.13) 
or (17.17) following a rotation of the type (12.10). It can be seen that 
this is precisely the sequence of transformations used by Chandrasekhar 
and Xanthopoulos (1987a) as described above. 

17.4 A technique of Halilsoy 

A simple technique was described in Section 16.2 in which a solution E Q 
of the vacuum Ernst equation (11.14) was used to generate a colliding 
electromagnetic wave solution of the equations (16.21) and (16.20) given 
by Z = 1, H = E Q . This technique was included in that of the previous 
section. In terms of the alternative Ernst potentials, the new solution 
satisfying (16.25) and (16.26) is given by E = 0, r\ = E Q . 

As pointed out by Halilsoy (1989a, b), this generation technique can 
easily be generalized. It can immediately be seen that a new colliding 
electromagnetic wave solution of the equations (16.25) and (16.26) can be 
given in terms of the alternative Ernst potentials by 



E = aE , rj = y/l -a 2 E Q (17.24) 

where a is an arbitrary parameter. 

It is also sometimes convenient to express this transformation in 
terms of an alternative parameter a such that 

a = cos 2a. (17.25) 

The new solution may then be expressed in the form 

E = cos2a£ , 7] = sm2aE . (17.26) 

In terms of the original Ernst potentials, the transformation (17.24) 
or (17.26) can also be stated in the equivalent form that, if Z Q is a solution 
of the vacuum Ernst equation (11.8), then a new electromagnetic solution 
satisfying (16.21) and (16.20) is given by 



Z 
H 



[1 — a) + (1 + a)Z Q sin a + cos 2 a Z Q 



(1 + a) + (1 - a)Z cos 2 a + sin 2 a Z Q (17 27) 

Vl — a 2 (Z — 1) sin a cos a (Z Q — 1) 
(1 + a) + (1 - a)Z Q cos 2 a + sin 2 a Z Q 



It is of interest to note the similarity between the initial forms of these 
expressions and those of (17.9). 
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Proceeding in a similar fashion as in (17.10), it can be seen that the 
potentials for the new solution are given in terms of the vacuum potentials 
Z = ^ + z$ by 

4 4a 
* = $ = 0^° (17 - 28) 



where, as in (17.9), 



2 = |(1 + a) + (1 - a)Z Q \ 2 
= 4| cos 2 a + sin 2 a Z \ 2 . 



(17.29) 



The expressions (17.27) may be substituted directly into (16.28), 
which then reduces exactly to (11.26). It follows that, if the initial solu- 
tion Z Q satisfies the boundary conditions for colliding plane gravitational 
waves, then the new solution (17.25) also satisfies the required boundary 
conditions for the same metric function U given by 

U = U Q . (17.30) 

It can be seen that the vacuum solution can be recovered when a = 1 
or a = 0, and that the original technique of Section 16.2 which is given 
by (16.33) or (16.34) is obtained when a = 0ora = 7r/4. 

The transformation (17.24) was first used by Halilsoy (1989a) to ob- 
tain a new class of solutions with the familiar form for the vacuum solution 
given by 

E Q =pt + iqz (17.31) 

where, as in (13.4), p 2 +q 2 = 1. The resulting solutions have been analysed 
in detail. In the limit as a — > 0, this solution must reduce to that of 
Chandrasekhar and Xanthopoulos (19866) as described in Section 13.3, 
which is the Kerr solution in the interaction region. 

It is also of interest to note that it is a very simple matter to generalize 
the transformation (17.24), or (17.26), that has been used in this section. 
It can easily be seen that, if E Q is a solution of the vacuum equation 
(11.14), then a family of solutions of the electrovac equations (16.25) and 
(16.26) is given by 

E = aE Q + b, n = cE + d (17.32) 

where a, 6, c and d are arbitrary complex constants that are only required 
to satisfy the two conditions 



ab + cd = 0, aa + bb + cc + dd = 1. 



(17.33) 
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It can immediately be seen that this generation technique is not only a 
generalization of (17.24) and hence also of (16.33), but is also a gener- 
alization of the transformation (17.8) that was described in Section 17.1 
and applied in Section 17.2. The generation technique given by (17.32) 
may also be expressed in the equivalent form 

Z=i±A±^, H= d t CK ■ (17.34) 
l-b-aE ' l-b-aE y 1 

By substituting (17.32) into (16.30), it can be seen that the boundary 
conditions for colliding plane waves are automatically satisfied for the new 
solution if they are satisfied for the initial solution E Q . It follows that, 
using the transformation (17.32), a large class of solutions describing col- 
liding combined electromagnetic and gravitational waves can be obtained 
from any colliding gravitational plane wave solution. 



17.5 Other solutions 

It has been pointed out by Neugebauer and Kramer (1969) that, for any 
space-time having a non-null Killing vector, there are a number of in- 
variance transformations for Einstein-Maxwell fields on the assumption 
that the electromagnetic field shares the space-time symmetry. These 
transformations can be considered to arise from the existence of scalar 
potentials and the symmetries of the associated potential space. This is 
well reviewed by Kramer et al. (1980). It is appropriate at the beginning 
of this section simply to summarize a number of results which follow from 
this, and which apply specifically to the colliding plane wave problem in 
which there is an additional space-like Killing vector. 

In terms of the Ernst potentials introduced above, the invariance 
transformations may be stated in the form that if (Z Q , H Q ) is a solution 
of the Ernst equations (16.22), then other solutions (Z,H) are given by 

Z = aaZ , H = aH Q , (17.35a) 

Z = Z + ib, H = H , (17.356) 

Z ° H= - H ° , (17.35c) 



1 + icZ Q 1 + icZ Q 

Z = Z Q + 2pH Q + H = H Q + (3, {VI Md) 

Z= ^ H= H ° + 1 Z ° ( i 7 . 35e) 

l + 2 7 tf + 77 Z ' 1 + 2 7 tf + 77 Z v ; 



where a, /3, and 7 are complex constants and b and c are real constants. 
Together, these constants form the eight real parameters of the isometry 
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group of the potential space. Clearly, these transformations can be used 
in any combination to generate new solutions of the Ernst equation. 

It may immediately be observed that, in the absence of the electro- 
magnetic field, (17.356) is the transformation (12.17) and (17.35c) is the 
Ehlers transformation (12.13). Also, if the initial solution is purely grav- 
itational so that H Q = 0, then any combination of these transformations 
which includes (17.35e) will generate electromagnetic solutions from the 
initial vacuum solutions. Applied in this sense with H Q = 0, (17.35e) is 
referred to as a Harrison transformation (Harrison, 1968). 

It can also be shown that, if the initial solution (Z ,H ) satisfies 
the boundary conditions (16.28) for colliding plane waves, then this con- 
dition is also satisfied by a new solution (Z, H) obtained by any of the 
invariance transformations (17.35) for the same values of the parameters 
ki. Thus, for any initial solution which is a colliding plane wave solution, 
further colliding plane wave solutions can be generated by these transfor- 
mations. This result is very remarkable since, in the context of stationary 
axisymmetric space-times, these transformations do not always preserve 
asymptotic flatness and so are not of great practical importance in that 
situation. 

It was initially pointed out by Kinnersley (1973) 2 that the group 
of symmetry transformations of the Einstein-Maxwell equations with a 
non-null Killing vector is the group 577(2,1). This is the group of the 
transformations (17.35). It has been further described in terms of colliding 
plane wave solutions by Li and Ernst (1989). 

It can now be seen that the approaches of Chandrasekhar and Xan- 
thopoulos (1985a) described in Section 17.1 and of Halilsoy (1989a) as 
described in Section 17.4 are both equivalent to 'Ehlers-Harrison' trans- 
formations for particular choices of parameters. 

Clearly a large class of further colliding plane wave solutions can be 
obtained using the above transformations. A four-parameter generaliza- 
tion of the Nutku-Halil solution has already been obtained by Garcia-Diaz 
(1988) using this method. Garcia-Diaz (1989) has also used a Harrison 
transformation to obtain a four-parameter generalization of the Ferrari- 
Ibanez solution. Further solutions have been generated in this way by Li 
and Ernst (1989) generalizing the solutions of Ernst, Garcia-Diaz and 
Hauser (1987a, 6, 1988), Ferrari, Ibanez and Bruni (1987a, b) and the 
Nutku-Halil solutuion. 

There is also a large number of more general solution-generating tech- 
niques that can be applied to this situation. These are effectively gener- 
alizations of those mentioned briefly in Section 12.6, which were initially 



2 See also Kramer et al. (1980) section 30.3. 
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developed in the context of stationary axisymmetric space-times. For ex- 
ample, Cosgrove (1981) has found a Backlund transformation, expressed 
in terms of the Hauser-Ernst formalism, which applies to electrovac space- 
times with two Killing vectors. 

It is finally appropriate here to refer to other non-diagonal solutions 
which describe colliding plane electromagnetic waves. 

A five-parameter family of type D electrovac solutions has been given 
by Papacostas and Xanthopoulos (1989). In these solutions, the hyper- 
surface / + g = in region IV is a quasi-regular singularity. The solution 
is in fact a special case of a family of solutions given by Debever (1971), 
Plebanski (1975) and Plebahski and Demiahski (1976). 3 There is also a 
vacuum limit of this solution which describes a collision of gravitational 
waves. 



See also Kramer et al. (1980) section 19.1. 



